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SUMMARY 


This report presents several examples in which flutter Involving 
Interaction between flight mechanics modes and eiastic wing bending occurs 

• I 

for a forward swept wing flight vehicle. These results show the basic 
mechanism by which the Instability occurs and form the basis for attempts 
to actively control such a vehicle. This work was accomplished between 
1 Nay 1981 and 31 October 1981. 
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INTRODUCTION 

% 

It has been known for many years that the calculated flutter speed 
of an aircraft may be modified by the Inclusion of the rigid body or 
free-free modes of the aircraft. Similarly, the predicted aircraft 
stability and control characteristics may be significantly different If 
structural flexibility Is Included In the analysis. Over the years, cer- 
tain Informal guidelines have been developed to assess the potential 
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Importance of the body-freedor*» modes to flutter analysis. 

The earliest report on the subject of body- freedom Influence on 
flutter appears to be due to Frazer and Duncan (Ref. 1). Their study 
showed that "fuselage mobility" had a slightly favorable effect on the 
calculated flutter speeds, Broadbent (Ref. 2) discussed the necessity 
for the Inclusion of the body-freedoms for swept wing flutter calculations. 
Gaukroger (Ref. 3, 4) also conducted several significant Investigations 
Into the Influence of body-freedom on flutter. His studies are summar- 
ized in the A6ARD Manual on Aeroelasticity (Ref. 5). His studies revealed 
that, when body-freedoms are Included In the analysis, two distinct types 
Of symmetric flutter can occur. The first type of flutter results from 
wing bending coupling with body pitch. Gaukroger uses the term "body- 
freedom flutter" to describe this phenomenon. He further notes that the 
flutter speed for this type of flutter may be low when compared to the 
flutter speed computed without body freedoms Included. In addition, the 
reduced frequency at which body freedom flutter occurs Is also relatively 
low. Measured In terms of wing chord lengths traveled by the aircraft 
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during one 05d nation cycle, this travel distance at the flutter speed 
is of the order of 100 chord lengths per cycle. These oscillations are 
well within the range of interest of flight mechanics specialists. 

The second type of symmetric flutter that can occur with body-freedom 

* f 

included is relatively insensitive to fuselage inertia. This type 

of flutter resembles that obtained when the wing root is held fixed. The 

$ * 

critical speeds, mode shapes, and frequencies are very nearly the same 
for the body- freedom and fixed root analysis. Gaukroger also determined 
that the presence of the horizontal ’Stabilizer had a favorable effect on 
the body-freedom flutter speed. Gaukroger* s results were confined to the 
aft sweep region. 

An additional, definitive, report on the subject of body-freedom 
flutter is given in Ref. 6 by Cunningham and Lundstrom. In Ref. 6, 
an analysis is presented to predict the flutter speed of a sounding 
rocket with an unswept, wing. Since their results have a direct bearing 
upon the results presented later in this report, we will review them here. 

The Cunningham/Lundstrom structural dynamic model of the sounding 
rocket included four degrees of freedom. 

(1) Wing bending in the fundamental mode. 

(2) Wing torsion in the first mode. 

(3) Fuselage pitch 

(4) Fuselage plunge 

Unsteady aerodynamic effects are introduced via the familiar Theodorsen 
aerodynamic theory used in conjunction with aerodynamic strip theory. 

Their analysis found two critical flutter roots. The first was a 
body- freedom root, while the second was a wing bending-torsion root. 

The former occurred at a frequency of 26% of the bending frequency, while 
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the latter occurred at 3.^40 times the fundamental bending frequency. A 
significant result of their analysis was that wing position relative to 
the missile c.g. had an important effect on the value of the flutter velo- 
city. As the wing leading edge approached the c.g. the body-freedom flut> 
ter speed was increased. -As the wing moved aft> the body-freedom flutter 
speed first declined, then slowly increased. 

t • 

In Ref. 6, the importance of including vari^>us degrees of freedom 
was discussed. From the calculations perfomed by modeling the system 
as a binary or ternary system, it was concluded that the inclusion of 
rigid body plunge was important, in that It raised the flutter speed when 
compared to a similar computation excluding this degree of freedom. In 
addition, although wing bending stiffness rather than wing torsional 
stiffness was found to have the predominant effect on body- freedom flutter, 
Including torsional flexibility lowered the flutter speed by about 14^. 

Since static aeroelastic divergence can be characterized as flutter 
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at zero frequency, we should be alert to the potential effects of body 
freedoms upon divergence. This problem Is lucidly discussed by Hancock 
(Ref. 7) for the case of torsional divergence; the effect of static aero- 
elastic flexibility upon the aircraft stability and control problem is 
discussed in Ref. 8, also by Hancock. However, until recently the forward 
sweep bending divergence problem has not been considered. 

The purpose of the present study is to attempt to identify the 
importance and the effects of the inclusion of body freedoms into a dyna- 
mic aeroelastic analysis of a symmetrical forward swept wing. For this 
study, we will employ a ternary model with body freedoms in pitch and 
plunge and a single fundamental bending mode. Torsion is deleted to iso- 
late the effect of bending flexibility on the forward swept wing divergence 
problen. Our focus will be confined to the low frequency body freedom 


flutter phenomena. The particular interest is the determination of where 
in the flight envelope body freedom flutter occurs relative to clamped 
wing divergence. For this reason, quasi -steady strip theory aerodynamics 
are employed since their mathematical representation is simple and rela- 
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tively easy to interpret. 

Model Parameters and Restrictions 

The planform geometry of the aeroelastic stability analysis model is 
illustrated in Figure 1. The flexible wing has a uniform chord, c, and 
sweepback. A. The junction between the wing and fuselage is located a 
distance (denoted as the dimension x) aft of the aircraft c.g. The swept 
semi-span dimension of each wing is ji. 

The model has three degrees of freedom; these are; 

1) vertical displacement of the aircraft c.g.. denoted as w, 
measured positive upward. 

2) pitching of the fuselage about the c.g., denoted as e, 
positive nose-up. 

3) bending of the wings, measured with respect to the wing root, 
positive upward. 

Motion is restricted to be synmetrical with respect to the fuselage center- 
line. This approach is essentially the same as that used by Pai and Sears 
(Ref. 9) to study the effect of wing bending upon aircraft longitudinal 
oscillation and by Cunningham and Lundstrom in the previously cited work. 

In the present analysis, however, the effects of an active canard/tail are 
included, since the presence of active controls can modify our conclusions. 

Strip theory aerodynamics are used to describe the airloads present 
on the wings as they oscillate. These loads are assumed to be quasi -steady 
so that the aerodynamic influence coefficients are independent of reduced 
frequency. 
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Three equations are necessary to describe the motion of the aircraft 
model, These equations are developed In the Appendix. One equation des- 
cribes the vertical motion of the aircraft center of gravity. A second 
equation describes the oscillatory behavior of the flexible wing, while 
the third equation Is written to describe the pitch motion of the air- 
craft. These three equations contain differential operators and Integrals 
(see the Appendix for details). By using Galerkin's Method, and suitable 
nondlmenslonallzatlon, the equations can be written In matrix form. Sev- 
eral parameters arise from this nondlmenslonallzatlon and will be discussed 
before proceeding further. 

Body-freedom Flutter Parameters 

With the assumptions and restrictions discussed previously, a set of 
parameters, some of them nondimenslonal, appear In the equations of motion. 
These parameters are as follows: 

u * 2 mif Mp “ the wing/fuselage mass ratio. Note that the 
total mass of the aircraft is Mp(l + y). 

V X + sInA = the longitudinal distance from the aircraft c.g. 

^ A 2 to a wing center span, measured In swept wing 
lengths. See Figure 1. 

Mp * Mp + 2 m£ = total mass of the aircraft. 

* (1 +y) Mp 
2 

Iq * Mp r * the aircraft mass moment of Inertia In pitch, 

^ measured about the aircraft c.g. 

• 2 r ^ Q 

■ aircraft radius of gyration, squared, divided by t . 

S * c£ - wing area of one wing 
Cj^a " 11^1^ curve slope 

w * natural frequency in bending for the wing with the 
fuselage fixed. In radians per second. 
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» >1 P cos^ A, dynatnic pressure parameter. 

V * V cos A. the air speed component normal to the vrtng 
" leading edge. 

With the assumption that the time dependence of the motion is the 
equation for the dynamic response of the model is ' 

[s*[H] ♦ s[B] + [K]] (t) - (F) (1) 

where 

C, * w/i 

C2 * V* 

?3 = 0 

and w» F. and e are modal amplitudes of the three distinct types of motion. 
The three matrices, [M], [B] and [K] are defined in the Appendix. 

The mass matrix [M] is nondimens ional and contains cross-coupling 
terms. These coupling terms are functions of the mass ratio p, the sweep 
angle A and the longitudinal distance between the aircraft c.g. and the 
single wing semi-span. One would suspect that, as in all flutter problems, 
the mass coupling plays an important role in either controlling flutter or 
triggering flutter. In particular, the coupling term M23 describes the 
inertia coupling between wing bending and fuselage pitch. It will be 
shown later that this term can be positive, negative or zero, depending 
upon the wing longitudinal position x and the sweep angle A. 

The aerodynamic damping matrix [B] is symmetrical. It is composed of 
nondimensional terms that are functions of x and sin A. It, too, is fully 
populated. The damping matrix has a common factor, called D in the analysi 



the stiffness matrix is composed of the structural stiffness of the 
wing with respect to the fuselage and aerodynamic stiffness tenns. The 
structural stiffness consists of a single term in the diagonal position 

associated with the wing bending mode deflection R. The aerodynamic 

• » ' 

stiffness matrix Is a function of sween, a. wing position x and a common 
factor* Q. 

# 

The aerodynamic stiffness matrix has several interesting features. 

The first column, associated with aircraft plunge, consists entirely of 
zeros. This Is consistent with the fact that a small change in altitude 
does not alter the lift. If only the Kg2 retained, corresponding 

with the assumption that the fuselage is fixed, then static aeroelastic 
divergence will occur when the wing is swept forward. This divergence 
occurs at 

, . rM56Er ,,, 

^La ^ ^ 

The above coefficient of -6.356 Is in excellent agreement with the exact 
solution for which the coefficient is -6.33. The difference occurs be- 
cause the assumed bending mode shape used in the analysis is a polynomial 
and not the exact solution. 

In addition, the fact that the aircraft system stiffness matrix has 
cross coupling terms proportional to dynamic pressure reminds us th.it the 
classical wing divergence problem must be modified in the case of the 
freely flying aircraft. Aeroelastic divergence will occur when the aero- 
dynamic loads cause the stiffness matrix to become singular. This subject 
will be discussed later, but two points will be made now. First, in our 
model, the plunge degree of freedom plays no role In the determination of 
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aircraft divergence because of the column of zeros associated with this 
degree of freedom, This degree of freedom is always neutrally stable 
with retivect to small perturbations as long as v remains constant. Second, 
the presence of the canard/tall alters the problem, as does wing placement > 

• X 

This is so because these parameters introduce cross>coupl1ng in the stiff- 
ness matrix and these terms appear in the characteristic equation used to 

* I ♦ 

determine divergence. 


The Model for an Active Canard 


A model of an all -movable canard was also developed. The geometry of 
this canard Is shown in Figure 1. The effect of this canard is assumed to 
be aerodynamic only. Any dynamic or inertial effects are ignored. With 
quasi -steady aerodynamic theory the canard aerodynamic forces and moments 
appear explicitly only in the c.g. and pitch equations of motion. The 
quasi-steady angle of attack of the canard surface is: 


o. ' e + «. 


- (^ + ^) 
' u • u * 


(5) 


It is recognized that the aerodynamic interaction between the canard 
and the wing is considerably more complex than that just described. How- 
ever, this approach provides at least a starting point for the discussion 
of active control of the aircraft with a canard. 

The presence of a canard/tail control surface will add terms to the 
right hand side of our equations of motion. These terms fall within two 
categories: terms dependent upon the aircraft degrees of freedom (c^); and 
terms depend upon the canard deflection, 5^^. The terms that are dependent 
upon the aircraft freedoms must be transferred to the left hand side of 
the equation of motion. Thus, the canard will now influence our aeroelas- 
tic stability problem, even when it is not active. The term proportional 
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to 6^ can be made to be a function of (c) by defining a feedback control 
law. This potential for improving the vehicle aeroelastic stability will 
be discussed later. 

I 

The important parameters for the canard are called 3 and f ;ln the 

• I 

analysis. These terms are defined as: 


3 ■ d/i (see Figure 1) 


f * 




where is the effective canard area and *^Loc 
slope. 


is the canard lift curve 


The Effect of Wing Position Relative to the Aircraft C.6. 


To Illustrate the potential effect of wing position Upon body freedom 
flutter, we will define an aircraft configuration typical of a fighter air- 
craft. It should be noted at the outset that the parameters chosen are 
typical of suggested designs for FSW fighters. Any resemblance to a 
specific design is coincidental. The parameters chosen are as follows: 


^Lo 

(1/rad. ) 


a 

f 

4M 

“o 

(rad/sec. ) 

i 

(ft.) 

My/2S 

(slugs/ft^) 

A 

(deg.) 

6.28 

0.11 

0.3 

0.17 

0.61 

68.0 

15 

3.8 

-30 


We begin our study by placing the wing in a far aft position with 
respect to the c.g., x « 0.45. We will display our eigenvalue analysis 
results in root-locus format for clarity. The “gain" in this case will be 
dynamic pressure, q. No attempt will be made to correct the aerodynamic 
representation for Mach nwiber effects, even though speeds are considerable. 

Let us first focus on the eigenvalue root that begins as the pitching 
or ''short period" mode. This root begins at the origin since, without 




I 


dynamic pressure, the pitching mode Is neutrally stable. A root-locus 
plot of this short period root Is shown In Figure 2 for a number of dyna- 
mic pressure values. The root 1r> oscillatory and well damped until a 
velocity equal to about 0.8 Vqq, where Is the clamped wing divergence 
speed. Up until this velocity, this mode has been acquiring a significant 
amount of bending deflection. After 0.8 Vqq, the root begins to move 

t * 

rapidly towards the right half plane as velocity is Increased, At V * 0.89 
the root becomes unstable at a frequency of 22.37 rad/sec. correspond- 
ing to a reduced frequency of k ■ wb/V * 0.026. 

The examination of the flutter mode reveals that at the time that the 

t 

bending displacement reaches its maximum upward value (R = 1.0), the air- 
craft plunges downward 0.55 units while the aircraft pitches upward 9.8 
degrees. Upward fuselage motion lags upward bending by 174 degrees while 
nose up pitch leads by about 8 degrees. 

Accounting for different coordinate definitions, these results are 
similar to those observed by Cunningham and Lundstrom for their rocket. 

Turning to the bending mode behavior, we see in Figure 3 that this mode 
Is more highly damped than the pitch mode. At about 0.8 Vp^ this mode 
begins to develop highly damped oscillatory characteristics. An increase 
in dynamic pressure is seen to drive the frequency of the bending 
mode down, while the pitch mode frequency Is increased. A classical fre- 
quency coalescence occurs, resulting in one mode becoming unstable. As 
velocity Is increased slightly beyond 1.0 Vpp, the bending mode contacts 
the real axis In Figure 3. After this point, the eigenvalues for this 
mode are no longer complex conjugates, but are real and distinct. One 
root moves left, in the direction of increasing stability, while the other 
rapidly moves right with increasing velocity. As this root enters the 
« right half plane at the origin, aircraft divergence occurs. 

i 

I 
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PITCH BRANCH 
X - 0.A5, A= -30" 


LABELLED POINTS 
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AIRCRAFT DIVERGENCE 
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Unlllice clamped wing divergence, where the aircraft is restrained 
against motion, aircraft aeroelastic divergence involves wing bending 
and fuselage pitching. Mathematically, the eigenvalue problem that deter^ 
mines aircraft aeroelastic divergence comes from the requirement that, at 
divergence, the total stiffness matrix, including aerodynamic and structural 
effects, must singular. The eigenvalues are related to the dynamic 
pressure. In terms of an eigenvalue solution for the dynamic equations, 
this divergence solution is valid if s = 0. 

For a particular planform geometry, and with given aerodynamic charac- 
teristics and structural stiffness (in our case, this stiffness is charac- 
terized by the value of the root-fixed bending frequency), the three 
degree of freedom model will yield three values of q for which the system 
is neffitrally stable. One of these values is zero, the second is unde- 

<r f 

fined, while the third value is a function of the previously mentioned 
variables.. 

The zero flight dynamic pressure value indicates that the aircraft 
is neutrally stable in pitch if there is no forward speed. This occurs 
because an aerodynamic moment about the c.g. provides a "pitch spring" 
that, in turn, provides the restoring moment if the system is perturbed in 
pitch. In addition, the system is neutrally stable in plunge at any dyna- 
mic pressure since there are no "plunge springs", aerodynamic or otherwise. 

The third root is meaningful. With the definition that q^Q is the root- 
fixed divergence speed, the aircraft aeroelastic divergence speed is q^^, 
given by: 


**DA 


q„ [ 5(y cos , A - if) 


] 


(y cos A -5 3f -0.4 sin A cos a) 
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M 

^DC 




S(y cosA - 3f) 


1/2 


(y cos A -5 3f - 0.4 sin A cos h] 


vihere y x 4 . 

There appear to be two Timitlng cases for the above formula,' If a 
tail is used* then 3f is less than zero. As 3f increases, due either to 
the use of a large horizontal stabilizer or a long tail length, the ratio 
Vda/Vdc approaches unity. On the other hand, when 3f is greater than zero, 
ratio of Vp^/Vp^, is infinite when 

5f k Z_|2i * 1^ sin A cos A 

As df increases from this value, the aircraft divergence dynamic pressure 
becon>es negative. At some value of 3f the airplane will be statically 
unstable. In this case, Vp^/VpQ is zero. As 3f increases, the ratio 
becomes more positive and finally approaches unity as df gets larger. 

These limiting trends are identical to those observed by Zeiler in 
a more general theoretical study of FSW aircraft divergence. It is clear 
that the analysis of static instability must reflect the presence 'of a 
canard/tail and include the aircraft geometry. 

The mode shape at aircraft divergence is given as: 

0 .. / V + TO sin A^ _ , . _ __ . 

— a — 1 Sin A cos A 

R y cos A - 3f 

If the aircraft is stable, then the numerator is greater than zero. The de 
nqminator will be positive if the airplane is statically stable with the 
canard in place. In this case, a FSW aircraft Will have an aircraft 
divergence mode consisting of upward wing bending and nose-down pitch 

(e/fi < 0). 
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For the aircraft whose dynamic stability properties were 
discussed previously, the ratio is 2,59, well above 

the clamped wing bending divergence speed. Thus, we can assume' 
that body freedom flutter is critical, unless a high frequency 
flutter mode is present at a lower speed. 

The Effect of Moving the Wing Root Forward 

Figure 4 is a root-locus plot of the results of an eigen- 
value analysis of the aircraft equations of motion. Both the 
bending mode and the pitching mode roots are illustrated as a 
function of velocity. In this example, x » 0.35 so that the 
wing root has been moved 0.1 OJi forward of that shown in the 
previous example. 

The striking difference between this example and that shown 
previously is that the "ancestry" of the unstable mode has its 
roots in the wing bending mode. Wing bending coupling introduced 
into the pitch mode stabilizes that mode and actually causes it 
to be highly damped. 

A similarity between the present case and the previous case 
is that, despite its ancestry, the flutter mode still occurs 
at about the same frequency and still involves coupling among the 
three degrees of freedom. Significantly, however, the flutter 
speed has increased by nearly 9 %. This case undoubtedly more 
closely represents an actual aircraft configuration with the wing 
center of pressure in close longitudinal proximity to the aircraft 
c.g. 
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Some parameter has apparently changed significantly to cause 

the flutter mode to switch its ancestry. The bending mode has 

been modified so that it no longer contacts the real axis in the 

left half plane. There is still an aircraft divergence speed, ' 

but it is at 6.08 and is associated with the pitching mode. 

One significant difference between the 7 ■ 0.45 case and 

the "x = 0.35 case is that the mass coupling term has changed 

sign, from minus to plus, as the wing root moved forward. The 

expression for this term in the mass matrix is 
» 

This term is zero if x “ where 

Xq * ^ sinA (11) 

For the 30 degree forward sweep case. 

Figure 5 shows the root locus plot for the case where x - 
The critical mode still appears related to pitch mode instability. 
However, a distinct "pinching together" of the bending root and 
pitching root is apparent. Comparing Figure 5 to Figures 3 and 4, 
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we see the transition of these characteristic root locus plots 
as the mass coupling changes. 

Figure 6 Illustrates the behavior of the bending root for 
)( ■ 0.30 as velocity changes. In this case^ the pitch mode Is ^ 

statically unstable and the roots for that mode lie 1n the right 

# 

half plane at all, values of velocity. Note that the flutter 
frequency Is relatively unchanged from other examples, but It has 
declined as x decreases. On the other hand, the flutter speed 
has Increased to a little over Z% above 

A summary of the figures presented previously Is contained 
in Figure 7. In this figure the ratio Vpp/VpQ Is shown as the upper 
curve and Is plotted on the left scale, against K while the ratio 
is plotted on the right, also against x, and is the lower 
curve. Two vertical lines are drawn to indicate: the position 
where mass coupling between bending and pitch changes sign; and, 
the position, x ~ below which the aircraft becomes statically 
unstable at any speed as a rigid aircraft. Mention should be made 
that static aeroelastic deformation causes the aircraft to be 
unstable at a low q before this value x. is reached. 

An alternative to changing the wing root position is to change 
the wing sweep at a constant root position. We would expect to 
see similar behavior of the flutter speed since the longitudinal 
C.P. position is changing. However, in addition, we are changing 
the aerodynamic coupling by changing the sweep. 

Figure 8 shows the results of a quasi -steady flutter study for 
our model where x « 0.40 and A Is varied from 0 to -50 degrees. 
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THE EFFECT ON FLUTTER SPEED AND 
FREQUENCY OF WING ROOT POSITION 
RELATIVE TO THE A/C C.G. 




THE EFFECT OF WING SHEEP FORWARD WITH 
if • O.AO, ALL OTHER PARAMETERS FIXED 


O 



SHEEP FORWARD ANSLE (DEGREES) 




In this figure, the ratio of the body-freedom flutter speed at 
a given sweep angle to the body-freedom flutter speed at zero 
sweep is plotted. This curve, Vp/VpQ» begins at 1.0 when A « 0 
and declines as the wing Is swept forward. The minimum value of 
* body freedom flutter speed is reached at A ■ -30 degrees. There- 
after, the speed 'increases slowly. 

As the wing is swept forward, the theoretical clamped wing 
bending divergence speed drops rapidly (it is Infinite at A * 0 
degrees) and reaches a minimum at A > -45 degrees. Also shown 
in Figure 8 is the ratio of Vp to Sc* This ratio begins at 0 
when A « 0 and increases with forward sweep angle. A vertical 
line marks the transition from rigid aircraft static stability 
to static instability. 

' SUMMARY 

We have seen that a forward swept wing aircraft can suffer 
from body freedom flutter due to coupling between bending and 
rigid body motions, in this case pitch and plunge. One signifi- 
cant parameter th determines the degree of severity of this 
type of flutter is the nondimensional distance X For far aft 
wing positions, mass coupling is introduced that tends to 
introduce wing inertia forces that cause nose up pitching motion 
as the wing vibrates upward. As the wing root is positioned closer 
to the aircraft c.g. this inertial coupling decreases and finally 
changes sign at a position "x^. If the wing longitudinal CP is 
close to the c.g. we will have an aeroelastic Instability occurring 
near the clamped wing divergence speed of the aircraft. However, 
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this instability is oscillatory, not aperiodic. While the reduced 
frequency of the body freedom flutter mode Is small, in the range 
of 0.01 to 0.03, the flutter frequency is of the order of one« 
quarter to one-third of the fundamental bending frequency of the' 
fixed root wing. 

How closely these results from a greatly simplified model 
match the actual behavior of a high speed aircraft with all its 
complexity remains to be seen. However, a simple wind tunnel 
test by Foist and Sanger at Purdue in April of 1981 showed that 
body freedom flutter of the type described could occur. In addition, 
although finite span aerodynamic effects and the effects of compres- 
sibility will surely change the values of the body freedom flutter 
speeds calculated, the ratio of flutter speeds computed to the 
fixed root divergence speed should be less affected. 

Finally, the significance of the fixed root wing divergence 
speed to the aircraft dynamic response characteristics should be 
discussed. 

Although the mode of instability is not fixed root wing divergence, 
the speed is highly meaningful and not to be ignored. This 
speed is a benchmark to which all wing flexibility effects must 
be referred. At about 0.6 V„c. the wing begins to rapidly lose its 
stiffness (aerodynamic plus structural). Beyond this point coupling 
between flexible modes and rigid body modes is likely to be 
significant. This coupling and decreased stiffness may have 
important consequences to handling qualities and fatigue life. 
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The system stiffness Is composed of the sum of the aerodynamic stiffness 
plus the structural stiffness. [fO Is computed to be: 



A Model for an Active Canard 


A model of an all -movable canard was also developed. The geometry of 
this canard is shown In Figure A. The effect of this canard Is assumed to 
be aerodynamic only. Any dynamic or Inertial effects are Ignored. The 
canard Induced aerodynamic forces and moments appear explicitly only In the 
c.g. and pitch equations of motion. The quasi-steady angle of attack of 
the canard surface Is: 



The canard lift per side Is: 

27 


( 4 ) 




This may be written as: ^ . 

•*" 4" S,c,^/scu. 

““ 5 a ^/jL Cr » ur/j[ 


The pitching moment, positive nose up. Is found, for each canard surface, 
to be: 


Me* 


( 6 ) 


The lift and pitching moment expressions are next nondimensionalized 
to conform to our nondimensional equations of motion. Now, contributes 
to the lift equation as follows: 
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> 

(8) 


(9) 

^'’ = If. 

The nondimenslonal pitching moinent contribution to the equations of motion 

Is: 

M a = L ^ ^3? J 

*■ 



where W ^ ^ | 

Cos*^jt ^ M-ri- 

(12) 

Bji « ^ Vo>*-^) 

Mr ^ 

(13) 

30 




The presence of the canard surfacp will add a set of matrix terms to the 
left hand side of the equations of motion. 



In addition, a term due to 6 appears on the right hand side. 



The canard deflection 6^ is related here to the three degrees of freedom 
in the problem by the following relationship: 

5c, - 07) 

where the terms C^. C2 and C3 are complex functions in general, for sim- 
plicity, let 


C,0') s + sCu -V ■s’-Cii, 

C.O') “ 

I 




The matrix equation of inotion now becomes; 

where M, § and ^ represent the augmented system matrices, given by; 



, • W - Ccf ' 

( 26 ) 



( 27 ) 


1 - [c-f 

Cci/i 

( 28 ) 


[B ] and [K ] include the canard damping and stiffness, respectively, with 
the control locked (6^ - 0)» 

The augmented mass, damping and stiffness matrices are: 



6 


SC| 


2q 




where D 

-$dCv-<53> 

■ "''n 5 V«T 

+ S'-f 6/a.A 

Finally, 

r1 

the augmented stiffness matrix is: 

f “OCno Q 

-S^0jo-$f 




“2. 






-l‘?/<« 


ofi -ii-jf 

I mA ^ . 
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